Finite size scaling and triviality of 4 theory on 

an antiperiodic torus 

Matthijs Hogervorst a ' fe and Ulli Wolff"'* 

a Institut fur Physik, Humboldt Universitat, Newtonstr. 15 
12489 Berlin, Germany 

b Departement de Physique, Ecole Normale Superieure, 24, rue Lhomond 

75005 Paris, France 



Abstract 

Worm methods to simulate the Ising model in the Aizenman random cur- 
rent representation including a low noise estimator for the connected four 
point function are extended to allow for antiperiodic boundary conditions. 
In this setup several finite size renormalization schemes are formulated and 
studied with regard to the triviality of cp theory in four dimensions. With 
antiperiodicity eliminating the zero momentum Fourier mode a closer agree- 
ment with perturbation theory is found compared to the periodic torus. 
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1 Introduction 



It is generally believed that the quantum field theory of self-coupled scalar fields 
in four dimensions is trivial. Then all effects of interaction terms go away in the 
true continuum limit. This is certainly true in perturbation theory due to the 
universal positivity of the perturbative /3-function at small coupling (see section 
I3.2l for more details). In this framework one also understands that there is a range 
of values for the mass and self-coupling such that in the effective theory at finite 
cutoff there can be both substantial interaction and only tiny cutoff effects in 
physical quantities referring to energies much smaller than the cutoff. This is the 
reason why the appearance of a scalar Higgs field in the standard model may not 
really be in conflict with triviality. Turned around, triviality even implies order of 
magnitude bounds for parameters in the Higgs sector. 

There also is the logical possibility that scalar theories could have sectors for 
which perturbation theory is simply irrelevant. In a series of papers that began 
with [TJ Liischer and Weisz have provided a lot of evidence that in the standard 
lattice formulation of the theory there is no such nonperturbative sector for low 
energy observables, which was also confirmed early on by Monte Carlo simulation 
[2]. In most numerical investigations the four dimensional Ising model is simulated 
which arises as the infinite bare coupling limit of 4 theory. While it is plausible 
that this is the 'least trivial' case the main reason for this choice is the availability 
of particularly efficient simulation methods like the cluster algorithm employed in 

m. 

Triggered by [3] a further boost in efficiency could be realized for Ising simu- 
lations. The simulation strategy of Prokof 'ev and Svistunov consists of sampling 
certain strong coupling graphs to arbitrary order instead of spin configurations. 
For arbitrary finite systems the expansion converges and the two representations 
are equivalent. The authors propose a simple update scheme for the graph en- 
sembles which does practically not suffer from critical slowing down. This was 
elaborated in [3] by showing that there are in addition dramatically improved esti- 
mators available for two-point correlations. A further step was made in [S] where 
it was noted that the simulated strong coupling form coincides with the random 
current representation constructed by Aizenman [6]. He has proved correlation 
identities that allowed him to rigorously prove triviality for the Ising model with 
more than four dimensions (excluding however D — 4). A decisive achievement 
was a combinatoric construction of an estimator for the connected four-point func- 
tion that analytically subtracts the disconnected part and then could be bounded 
sharply enough. In the numerical approach at tracing the renormalized interaction 
strength [2] the corresponding cancellation had to be performed numerically which 
leads to much reduced precision. In [5] Aizenman's identity could be integrated 
into a numerical framework by simulating two replica and identifying the required 
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percolation clusters. It was demonstrated to lead to precise estimates of renormal- 
ized couplings with only moderate computing power in dimensions D = 3,4, 5. 

The question of triviality is an issue of ultraviolet renormalization. In numeri- 
cal simulations at D = 4 one is limited to ratios L/a < O (100) at present, where 
a stands for the lattice spacing and L is the system size. To investigate ultraviolet 
behavior close to the continuum limit a clever use of resources is achieved if one 
uses L itself as the scale to formulate renormalization conditions, i.e. to use a 
finite volume scheme as one does with the Schrodinger functional in QCD. Such a 
strategy - on a simple periodic torus - was in fact followed in [5] by fixing z = mL 
where m is the renormalized mass in the so-called second moment definition. It 
turned out that the first series of experiments at z = 2 led to a borderline agree- 
ment with perturbation theory: The leading order (one loop) described well the 
cutoff evolution of the renormalized coupling. The inclusion of two further orders 
of the asymptotic expansion however led away from the data. This situation was 
further investigated in [7] with the following results: At z = 4, closer to the ther- 
modynamic limit, perturbation theory works as expected yielding a successively 
improving precise description of the numerical data for the available three orders. 
For smaller systems it was argued that the single constant zero momentum mode is 
responsible for the bad 'convergence' of perturbation theory. This was supported 
by finding that an improved expansion, where this mode was treated exactly, yields 
a much better agreement with the data also at z = 2 and even for z — 1. 

In the present paper we have eliminated the constant mode by introducing an- 
tiperiodic boundary conditions in one or several directions of the four dimensional 
torus. Our new results here consist of generalizing the strong coupling/random 
current formulation to these cases including a proof of Aizenman's identities in 
section 2. In section 3 we define several finite size renormalization schemes and 
report values of the corresponding three loop /3-function coefficients. In section 4 
numerical results for the antiperiodic case are reported followed by some conclu- 
sions. Some details of the proof of the Aizenman formula and of the perturbative 
calculations are deferred to appendices. 

2 Random current form of the Ising model with 
antiperiodic boundary conditions 

2.1 Partition functions with charge insertions 

We here generalize the content of [5] while we at the same time slightly change the 
notation. We now start from the partition function 

Z[q\ = 2- L4 J2 e2Kl2l={xv)Z{l)s(xMv) n(s(x)) q(x) - (!) 
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We independently sum over s(x) = ±1 on all sites of a dimensional hyper- 
torus. Here and below all sites are understood to have integer coordinates with 
^ x M < L in all directions \l — 0, 1, 2, 3. On the nearest neighbor links I a Z(2) 
background gauge fielc{§ z{l) = ±1 enters and q(x) are fixed integer local charges 
whose values enter only modulo 2. For a collection of sites x^\x^ 2 \ . . . ,x^ n ' we 
may define 

n 

qi2...n(x) = ^2 $ x>x (i) (mod 2) (2) 
i=i 

and then a two point function is for example given by 

(s(x^)s(x^)) = (3) 

The sole reason for having the background gauge field here is to allow for antiperi- 
odic boundary conditions. We shall set 

= {xy)) = | (_"!)- ^ fe,V,} = {0,L-l} (4) 

Here e M is a 4-vector with zeros for the periodic and ones for the antiperiodic 
directions. Thus there is a minus sign on those links that 'close around the torus' 
in antiperiodic directions. This is equivalent to (partially) antiperiodic boundary 
conditions for the spin field. 

We next expand (CQ) in k by introducing an integer link field k(l) summed over 
values 0, 1,2, ...,oo (independently on each link) and, after summing over the 
original s(x), we arrive at 

Z[q] = y %2w[k]$ e [k]6 dk , q (5) 

k 

with 

»w=n%- w 

Here the divergence of k 

dk{x) = MO (mod 2) (7) 

l,dl3x 

is a site field where we add the k{l) of the eight links surrounding x. It is locally 
constrained to equal (modulo 2) the source q(x). The sign 



<m*] = n 



.lW)]* <0 - 
I 



1 We discuss D = 4 here, but the generalization to other D is trivial. 

2 We prefer z(l) over the usual notation z(x, ^) here because links are unoriented. 
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is the Z(2) winding number of the k field with respect to the antiperiodic directions 
coded into s. 



2.2 Connected four point function 

If we introduc^l 

Z c (<?i234) = Z[q 12U }Z[0] - Z[q 12 ]Z[q u ] - Z[q 13 }Z[q 24 ] - Z[q 14 }Z[q 23 ] (9) 
then the connected four point function is given by 

(s(x^)s(x^)s(x^)s(x^ c = (10) 

With the help of the results of appendix |A] this can be written as 

ZMizaa) = -2Y, w [ k } w [ k '}^M^e[k'}8 dktqi2 8 ak ', qM X lz , (11) 

k,k' 

where we have inserted 

<S> £ [k + k'] = $ £ [k}$e[k'} (12) 

for the arbitrary function F[k + k'). The constraint X 13 = X(x^\ x^; k + k') E 
{0, 1} is one if and only if x^ 1 ' and x^ are in the same percolation cluster with 
respect to bonds which are active on links where k + k! does not vanish. 
We now introduce an ensemble of two independent (factorizing) replica 

Z = Yl Yl W i k } W i k '] S 9k,q u Jdk',q u/vl (13) 
u,v,k u',v',k' 

with 

quv{x) = 6 X , U + 8 XiV (14) 

and expectation values ((•••)) in this ensemble are defined in the obvious way. In 
particular, the two point function now reads 

«->*»» = ff&fcff) (15) 

with the sign being part of the observables. This may of course be symmetrized 
over the two replica. For the connected four point function we get 

, ( / (4)nx o ( ( $ £ [k] $ £ [k'}5 xW i Jjy iV 5 x{3 ) y S xW V ,X (u, u'; k + k'))) 

{S[ x )-.-s { x )}e - L -8 (($ e [k]$ e [k']8 UtV 5 u > }V >)) 

(16) 

Note that for nonzero contributions, all four charges at u, v, u', v' are in the same 
percolation cluster. 



5 There is a trivial error in eq. (9) of [S]. 
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3 Renormalized mass and coupling 



3.1 Families of finite volume renormalization schemes 

Employing (partially) antiperiodic boundary conditions we now define a choice of 
possible renormalization conditions that use the finite system size as renormaliza- 
tion scale. The sets B £ of admissible momenta depend on e and are given by 

<t>(p) = e~ ipx s(x), P eB £ = { PlM = (n„ + ej2) x 2vr/L, < < L} . (17) 

X 

The zero momentum mode, which has led to a bad convergence of renormalized 
perturbation theory [7] on small tori, does not occur anymore for e ^ (0, 0, 0, 0). 

We use an index 'p' for the periodic case (eff 1 = 0), 'A' for fully antiperiodic 
(sjit = 1) and 'a' for one antiperiodic direction {s^ = 5 Mi0 ) an d then also write 
E s = B £ (s). For each case s G {p, a, A} we single out two small admissible momenta 
p s and p' s with distinct p 2 s < p' 2 with 

p 2 = 4^sin 2 (^/2), (18) 
n 

namely 

Ps = p s) , K=p s + (0,0,0,2tt/L). (19) 
Renormalized masses m s are defined by solving for m s in the universal ratios 

o _ (\m)\ 2 ) _ P 2 s+m 2 s _ 

Ks ~ 7TI7 — \m ~ ~ Z ?T\ 2 ^ z s — m s L. [20 

{\<f>{Ps)n Ps + mj 

We note that in the symmetric phase massive scaling regions defined by < m 2 <C 
1 we must adjust R p > 0, R a > 1/5 and R A > 1/3 (up to 0(L~ 2 )). 

To define corresponding renormalized coupling constants g s we employ the 
smaller of the two momenta and form another ratio 

Note the connected four point function in the numerator. The coincides 
with the scheme studied in [5]. Each definition of g s has the following properties: 

• It derives from universal renormalized ratios of correlations with wavefunc- 
tion renormalization factors canceling. 

• If applied to </> 4 theory away from the Ising limit, it coincides with the stan- 
dard bare coupling at tree level of perturbation theory. 
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• For z s — > oo boundary conditions become irrelevant and all g s coincide with 
the usual coupling defined by vertex functions at zero momentum, and m s 
approaches the infinite volume mass defined at zero momentum. We hence 
make contact with the scheme of pp. 

We may now parameterize the renormalized theory by fixing a mass and a coupling 
constant. In the most general case we may even choose different boundary condi- 
tions s,t G {p, a, A} for this purpose and approach the finite volume continuum 
limit L = L/a — > oo for fixed z t and g s . 



3.2 Beta functions 

For each set of normalization conditions we consider the Callan-Symanzik evolution 
equation with the cutoff L = L/a 



L^f\ Zt=z =-f3 sAz (g s ). (22) 



As is customary in the 4 literature pQ we take the derivative on the left hand side 
at fixed bare coupling. Terms of order L~ 2 are neglected and therefore (3 s ,t,z is a 
function of g s only. In the perturbative expansion 

/Ws) = E fo ?U m ( 23 ) 

the first two coefficients are scheme independent, 

b {1) b {2) -1111 (24) 

(47r)2' (4tt) 4 ' 1 } 

while the third coefficient is known in the infinite volume [TJ 

€L - (25) 



Knowing the two loop relation between couplings in two schemes allows to 
compute the difference between their respective coefficients (3 S / Z . All necessary 
formulas are found in section 3.2 of [7]. In the same paper, from Table 1, values 
for 6p 3 p, 2 (the scheme used there) for many different z are given by relating them 
in steps to large, effectively infinite, z. We use these values now and connect to 
them the schemes relevant in this paper. By working out the necessary Feynman 
diagram sums (see appendix IB. II for details) up to L = 100 we have obtained 
Table HJ 
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s, t, z 


/7 (3) 7 (3) \ / A \ A 


KL x ( 47r ) 


p,a,2 


237.2805 


646.3126 


a, a, 2 


-374.8514 


34.1807 


A, a, 2 


-377.5345 


31.4976 


p, a, 3 


4.65458 


40.9003 


a, a, 3 


-7.07988 


29.1658 


A, a, 3 


-6.20041 


30.0453 



Table 1: Three- loop coefficients of the /3-functions for the boundary conditions 
relevant in this study, i.e. mass z a = m a L = 2,3 and couplings g p , g a , qa- 



3.3 Estimators 

The correlations entering into m s and g s are now translated into expectation values 
in the ensemble (I13p . It is not difficult to find for the ratio in (|2"0"|) 

R s = iipk F7TT7 TTT) = * E («) (26) 

(($s[k]f a {u - v))) 

with 

fs(x) = Y[cos{Xf,p s ^), f s {x) = Y[cos{x^p' SII ), (27) 

H fj, 

The invariance with respect to separate reflections of each direction has been used 
to factorize the Fourier exponentials into cos factors^. Where e^ a \p a) p' a and p' A 
treat the four directions differently we average over all possible ways of putting 
the anisotropics in our observables. In addition we combine the two replica in the 
error analysis as discussed in [S]. 

For the coupling, a possible estimator is given by 

g s = 2(z 2 s +L 2 p 2 s ) 2 X s (28) 

with 

= ((<5> s [k}<5> s jk>}f s (u + u'-v- v')X(u, u'; k + k'))) 

((® s [kmk>]f s (u-v)f s (u>-v>))) 1 ] 



This is however a special choice. Since the left hand side of ( 116]) is symmetric in its 
arguments, one could also permute the arguments in f s (u+u'—v—v'), inequivalent 
choices being f s (u — u' + v — v') and f s (u — u' — v + v') in addition. We found 
it quite profitable in terms of errors to average over the three possibilities (after 



4 Omitted parts with sin factors would average to zero but still contribute noise. In principle, 
if it is anticorrelated with the signal, this could lower the error, but this is unlikely. 
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verifying that their mean values are compatible). The denominator factorizes of 
course in the replica, and also here other choices are possible^]. 

We note that the factors L 2 p 2 A f=s 4:L 2 p 2 ~ 4tt 2 will enhance the values of the 
actually measured observable quite significantly for the antiperiodic cases. Corre- 
spondingly X a A will typically be found much smaller than X p due to cancellations. 
The same cancellations - absent in X p which has a non-negative estimator - will 
also lead to lower achievable precision for the antiperiodic cases, in particular for 
i A\ The fluctuating sign included in our observables is however not of the kind 
that leads to an exponential signal to noise problem on large lattices. The signs 
are 'coherently' related to the winding around the torus and to Fourier modes 
with wave numbers of order 1/L and do not combine nearly independent signs 
from many small subvolumes. The latter is typical when 'sign problems' render 
numerical estimates impossible. 

3.4 Couplings from partition function ratios 

In our finite volume simulations we also obtain information about ratios of partition 
functions with differing boundary conditions, for example 

Za_ = (($A,tt)) /o m 

Z A «<M^)>' 1 ] 

From the form of these estimators where signs are averaged it is trivial that Z a ^ Z v 
and Za ^ Z v holds (in the Ising limit!), and in addition Za ^ Z a is plausible and 
indeed found numerically in all cases. In the Gaussian limit discussed below the 
same ordering holds. 

In perturbation theory we may expand the differences in free energy 

HZ t /Z s ) = f Q ' s (z t , L/a) + fl' s (z t , L/a)g + O(g 2 ) (31) 

and some values of f , fi are listed in appendix iBl 

We may hence define further conventionally normalized coupling constants by 

\n(Z t /Z 3 ) - fo' s 

h t/s = TtTs • ( 62 ) 

Ji 

They measure the response of the free energy to a change of boundary conditions 
and are expected to be physical quantities in the continuum limit and thus le- 
gitimate renormalized couplings in a finite volume. Note that the renormalized 
mass Zt is chosen on the right hand side of ( 13~T|) . To the order cited we may also 
just replace go by any renormalized coupling. Thus f^'l refer to relations between 



'The one given looks most natural and promising. 
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renormalized quantities and are expected and indeed found to reach finite limits 
as L/a — > oo. A numerical disadvantage of h t / s is the required subtraction of the 
tree level part. After this cancellation, the precision of h a /A in our present simula- 
tion is too low for a meaningful study of its evolution. Phrased differently, in our 
simulations Z a /Z^, which is typically per mil accurate, is given within errors by 
its free field value, another result that is consistent with triviality. 

4 Numerical results 

We have developed a serial C-code to sample the graphs of the ensemble con- 
tributing in f )13p . Details are very similar to those given in [5] except that we have 
this time used a Metropolis rather than a heatbath step to move the worm-ends. 
An iteration consists of L 4 (attempted) worm moves for each of the two replica 
with about 64 percolation processes interspersed to compute X. The signs $ s are 
'updated' as the worms move and are hence available at any time to continuously 
accumulate observables during the updates. An iteration requires a computational 
effort proportional to L 4 that is roughly comparable to a sweep in a local update 
scheme. For each data point (value of L and k) we executed 10 6 iterations (after 
equilibration) where we always stored the blocked measurements from 10 succes- 
sive iterations. Thus we had to analyze time-series of length 10 5 and in these 
units we found at most integrated autocorrelation times of unity and for many 
observables the absence of any relevant correlations. Equilibration under these 
circumstances has been unproblematic. 

We have exploited some trivial parallelization to generate our data by running 
between 4 and 64 copies of the system to produce the total statistics. The runs 
took place on dual-quad-core X2270 PCs. Each of the L = 64 runs took about 
3000 core-hours or about 2 days with 64 cores and the smaller lattices follow by 
scaling proportional to L 4 . 



L 


2k 


Z a 


ft Clz a 
L 2 dK 


Zp 


z a 


za 


8 


.1475570 


1.9988(68) 


1.0804(52) 


2.2587(36) 


2.0 


1.781 (24) 


10 


.1482830 


2.0003(66) 


1.0440(50) 


2.2404(35) 


2.0 


1.835(21) 


12 


.1486864 


1.9948(64) 


1.0272(49) 


2.2270(34) 


2.0 


1.843(20) 


16 


.1490990 


2.0012(60) 


0.9739(47) 


2.2077(33) 


2.0 


1.843(18) 


22 


.1493687 


1.9971(57) 


0.9400(45) 


2.1947(31) 


2.0 


1.886(16) 


32 


.1495330 


2.0015(54) 


0.9043(42) 


2.1806(29) 


2.0 


1.884(14) 


64 


.1496509 


1.9993(48) 


0.8470(39) 


2.1611(26) 


2.0 


1.912(12) 



Table 2: Numerical results for simulations at z„ = 2. 



10 



9 P 


9a 


9a 




K /A 


16.44(7) 


30.26 (27) 


33 (12) 


16.58(52) 


33.9(6.6) 


15.14(6) 


27.19(25) 


36(11) 


14.25(50) 


23.1(6.3) 


14.25(6) 


24.94(24) 


48(11) 


14.10(48) 


31.8(6.1) 


12.96(6) 


22.45(22) 


19(10) 


13.26(45) 


19.7(5.6) 


11.89(5) 


20.23(19) 


7(9) 


12.10(42) 


21.9(5.1) 


10.87(5) 


17.81(17) 


36(8) 


10.99(38) 


16.9(4.6) 


9.34(4) 


14.70(13) 


11(8) 


10.09(31) 


11.5(3.7) 



Table 3: Companion to Table [2] with more observables, lines in the same order. 



In Tables [2] and [3] we compile results where for a series of lattices of sizes 
L = 8, . . . , 64 we have tuned k to values that produce z a pa 2 to a very good 
approximation. For all our observables we have measured their ^-derivatives as 
connected correlation with = J2i k (0j f° r example 

K Tk " w )» = ^Wsiu - v)S k )) - (($ s [k]f s (u - v)))((S k )) . (33) 

The fourth column of Table [21 relevant for the tuning of k has for example been 
obtained in this way. Moreover we have implemented a small post-run reweighting 
to first order in the K-shift to achieve z a = 2 exactly, as already discussed in 
[5]. The overall error estimate for this somewhat involved function of primary 
observables was determined following [8] with the error of the ^-derivatives safely 
neglected for the only small corrections. Thus the first four columns in Table [2] 
refer to the parameters that were actually simulated. The remaining columns as 
well as the couplings in Table [3] include the (tiny) corrections and thus refer to 
Zq = 2 as required for finite size scaling. We see relative errors fall with growing L. 
As we spend about constant computer time per site we experience slightly negative 
critical slowing down here. 

While the estimator ( 12 9p is non-negative in the periodic case this is not anymore 
the case for s = a, A. In particular for the fully antiperiodic coupling g& at z — 2 
the sign fluctuations are too strong to leave a useful signal at our statistics. 

Tables H] and [5] are structured in the same way as Tables |2] and [3] but refer to 
physically larger volumes with z a = 3. 

Remember that the columns in our tables refer to fixed values of z a and the bare 
4 coupling which is infinite in the Ising limit. For small a/L we expect (almost) 
universal relations between one pair (zt,g s ) and another one. The values of z p 
or za following downwards the columns of Table [2] or H] are expected to converge 
slowly at a rate given by the vanishing renormalized coupling (if triviality holds) 
to the free field result z p = z& = z a and not to nontrivial values at rates L~ 2 . 
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L 


2k 


Z a 


K C)z a 

L 2 dn 


Zp 


z a 


za 


8 


.1450850 


2.9968(42) 


0.7938(26) 


3.0815(24) 


3.0 


2.913(8) 


10 


.1465910 


3.0016(39) 


0.7593(25) 


3.0739(22) 


3.0 


2.919(8) 


12 


.1474720 


2.9980(37) 


0.7319(24) 


3.0661(21) 


3.0 


2.922(7) 


16 


.1483860 


2.9983(35) 


0.6946(23) 


3.0615(20) 


3.0 


2.942(6) 


22 


.1489732 


2.9969(34) 


0.6685(23) 


3.0550(18) 


3.0 


2.944(6) 


32 


.1493373 


2.9971(31) 


0.6337(22) 


3.0478(17) 


3.0 


2.962(5) 


64 


.1495982 


2.9988(28) 


0.5898(21) 


3.0417(14) 


3.0 


2.963(4) 



Table 4: As Table [2] but for z a = 3. 



9 P 


9a 


9a 


ka/p 


h a /A 


28.55(10) 


36.86(16) 


42.6(2.5) 


28.8(3.5) 


23(12) 


25.73(9) 


32.71(14) 


40.7(2.2) 


25.3(3.2) 


22(11) 


23.84(8) 


29.81(13) 


34.6(2.0) 


26.9(3.0) 


31(11) 


21.30(7) 


26.38(11) 


28.6(1.8) 


26.4(2.8) 


7(10) 


19.13(6) 


23.04(9) 


26.3(1.6) 


15.1(2.6) 


19(9) 


17.00(5) 


20.30(8) 


22.9(1.4) 


17.1(2.3) 


20(8) 


14.15(4) 


16.49(6) 


18.9(1.2) 


16.2(1.8) 


8(7) 



Table 5: Extension of Table HI 



In Fig. [U we see the evolutions of g a at z a = 2, 3. In both cases we find 
agreement with the perturbative pattern that systematically improves with the 
loop order for the three terms that are available. This is particularly pronounced 
if one compares the curves for z a = 2 with the fully periodic case in [5]. The 
expectation that antiperiodic boundary conditions render all modes perturbative 
also in a smaller volume is confirmed. 

In Fig. |2] we investigate two schemes in terms the periodic g p combined with 
z a . In particular the left plot looks very similar to [5]: the antiperiodic mass alone 
is not sufficient to eliminate the effects of the constant mode that contributes to 
9 V - 

In Fig. |3] we study two couplings which we can compute with lower but still 
reasonably significant precision. On the left we see agreement between Z a /Z p at 
z a = 2 with two loop perturbation theory, where no three loop term is available. 
Finally the right plot shows the fully antiperiodic g& at z a = 3. Since relatively 
large momenta contribute here, the cutoff effects could be larger in this plot than in 
the previous ones. While we expect them to be still small for L = 64, 32 we cannot 
really disentangle them. We refer to the discussion in [7] that the perturbative 
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2 loop 



1 0.2 0.25 0.3 0.35 0.4 0.45 0.5 1 §.2 0.25 0.3 0.35 0.4 0.45 0.5 

1/ln(L/a) 1/ln(L/a) 

Figure 1: Cutoff dependence of the coupling g a at z a = 2 (left plot) and z a = 3 
(right plot). 



artefacts, that could be specified at 1 and 2 loop order, are not relevant here. 

5 Conclusions 

In [5] a Monte Carlo algorithm was presented for simulating Aizenman's reformu- 
lation of the Ising model as a statistical system of random currents on links [6]. 
This form has two advantages: (practical) absence of critical slowing down and the 
availability of a non-negative estimator for the connected four point function, that 
enters into the standard definition of renormalized interaction strength, without 
having to perform numerical cancellations. In the present paper we have gener- 
alized this from periodic boundary conditions in all four directions to arbitrary 
combinations of periodic and antiperiodic directions. Then the above mentioned 
estimator starts to fluctuate in sign but still yields good precision for one an- 
tiperiodic direction and system sizes as small as z = 2. With four antiperiodic 
directions the noise was found to much more degrade the possible precision. One 
and the same simulation produces results for all the boundary conditions consid- 
ered here. For large z the simulated graphs do not wind around the torus and the 
independence of results on the choice of boundary conditions becomes manifest. 

In our previous work for the periodic case we have found, that the decay of 
the coupling strength with the UV cutoff (triviality) is not very well described by 
perturbation theory given the precision of our new method. This is a problem, 
because even very efficient numerical simulations cannot trace this decay over sig- 
nificant scale-ranges, if it happens only at the expected logarithmic rate. What 
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1/ln(L/a) 1/ln(L/a) 

Figure 2: Cutoff dependence of the coupling g p at z a = 2 (left plot) and z a = 3 
(right plot). 



it can only do is to confirm the matching with perturbation theory which can 
then be trusted all the way to the continuum. In [7] we have accumulated some 
evidence that the constant zero momentum mode that exists for periodic bound- 
ary conditions is the main source of nonperturbative behavior in small volumes. 
With at least one antiperiodic direction the smallest momentum is n/L and all 
modes receive Gaussian damping independently of the mass term. In Fig. [1] we 
demonstrate that perturbation theory indeed works much better now. 

In the present simulations we could determine with good precision the change 
in free energy caused by differing boundary conditions. Also these quantities can 
be related to renormalized couplings and probed numerically. As a nonvanishing 
tree level term has to be subtracted however, the precision is limited here. In per- 
turbation theory only the first two universal terms are known for these couplings. 
Given these limitations also here reasonable agreement is found (Fig. [3j left plot). 
Acknowledgements: We thank Peter Weisz for a careful reading of the manuscript 
This work has emerged from the internship project of M. H. at Humboldt Uni- 
versity and we thank the Ecole Normale Superieure (Paris) for financial support. 
UW. acknowledges support by the DFG via SFB transregio 9. 
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Figure 3: Cutoff dependence of the coupling h a / p at z a = 2 (left plot) and for the 
coupling gA at z a = 3 (right plot). 



A Proof of eq. ( HIT ) 



We freeze k (l) + k' (I) = K (I) to fixed values and then show for an arbitrary charge 
distribution p(x) the counting identity 

(34) 

where the large brackets are binomial coefficients. The sums here run over values 
from to K(l) independently for each k(l). If x and y are not in the same 
percolation cluster made from bonds with K(l) > 0, then X vanishes and the 
same is true for o~dk,q xy for all k ^ K, and thus both sides are zero. We thus 
only have to consider the case X — 1. Following Aizenman [B] (Lemma 3.2) we 
momentarily think of a graph C of K(l) distinguishable lines drawn 'over' each link 
with (link- wise) cardinality \C\ = K. Then on both sides we count the number 
of distinct subsets £' C £ which are such that the cardinalities k = \£'\ and 
K — k = \C\C\ obey certain divergence constraints. The identity is proven now 
by constructing a one-to-one mapping between the respective subsets contributing 
on the right hand side and the left hand side. 

Because of X = 1 there is some chain of specific lines C xy connecting x and y. 
If we now have some £' that contributes to the left hand side, then we map 

£' -)• £" = £'A£ xy = £'U £ xy \{£' n £ xy }. (35) 

One may now show that 
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if \C'\ satisfies the constraints of the left hand side of (151)1 then \C"\ fulfills 
those on the right hand side, 

under the same mapping C" maps back to £'. 



Therefore ( I34p is now established. 

We now multiply both sides of (134]) with where F is arbitrary and 

w is defined in (J6]) . Then we sum over K and change variables to k! = K — k and 
k to derive 

^F[k+k']w[k}w[k'}5 dKp 5 dk , Axy = y^Xfa y; k+k')F[k+k']w[k}w[k'}5 ak)P+qxy 5 dk > fl . 

b b 1 b b> 

i\j ^ r\i fx, . r\i 

(36) 

This may be used repeatedly to derive 



}] F[k+k')w[k)w[k') {5d k , qi234 5 dk/ ,o - S d k, qi2 S dk ',q 34 - 5 dk . qw 5 dk : m4 - 6dk, gi4 5dk',q 23 } = 

k,k' 

F[k + k']w[k]w[k'}S dkm234 5 dk ', o{l ~ Xm ~ *24 - X 23 } := E (37) 

k,k> 

with the short hand notation X^ = X{x^\x^\ k + k'). To satisfy dk = (/1234 and 
dk' = 0, the four points must either belong to one or to two percolation clusters 
of k(l) + k'(l) > 0. In the latter case it is easy to verify that {1 — X u — X 2i — X 23 } 
vanishes while in the first case it equals —2. This may now be changed back to 

E = -2 F l k + k'}w[k}w[k'}S 9km2 S dk 'm,X 13 . (38) 

k,k' 

B Perturbative expansion 
B.l Couplings based on correlations 

We here report details on the computation of the expansion (in go) coefficient^! 
Pi(s,t, z, L) and p 2 (s,t, z, L) for the schemes using g s , z t with s,t £ {p, a, A} de- 
fined in section 13.21 In terms of 1PI vertex functions our renormalized parameters 
read 

Y _, r^(p s ,-p s )-r?>(p> s ,-p> s ) 



( 2 )^ _„ ^ _ r( 2 ),„ _„ 



Ps - V 

,2 _ vr( 2 ) 



-Y t Y^\p u -p t )-pl (40) 
g s = -Y s 2 T^(p s ,p s ,-p s ,-p s ), (41) 



3 We count on our readers to not confuse these p with momenta. 
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with Y s being the 'wave function' renormalization factor and all quantities are at 
finite L. Standard bare perturbation theory gives (dropping subscripts s,t here) 

T (2) (P, -P) = -m 2 -f- g Q \ji + ^o^i^i(O) + 9l\Up) + O(g 3 ), (42) 

T^(p,p, -p, -p) = - go+ g^ Hl (p)-g 3 \3H 2>1 (p) + ~H 2 M + \h 2 ^p)\ +O(g 4 ). 

(43) 

The capital letters stand for the usual Feynman diagrams for the two and four point 
functions up to two loops and are given explicitly below. The mass parameter in 
all propagators is the bare mass wiq at this stage. The actual evaluation proceeds 
via the following sequence of steps, 

G(p)= * 2 , JV = 2simjv/2), (44) 



Note that in 



^)4E efp *> (45) 

qeB s 

Ju = G s (0). (46) 

G?(p) = Y.{ G ^ ne ~ iPX > ^ = 2,3,... (47) 



only integer momenta (e.g. 2p) are appropriate for even n where [G s (a;)] n is peri- 
odic for all s. We form (dropping subscripts s again) 

H 1 (p)= 1 -[2G^0) + G\2p)} (48) 

and similarly 

J 2 {p) = G3(p), (49) 

#2,i(p) = ^lE G ^)G 2 {q - p)[2G(q) + G(q - 2p)], (50) 

qeBs 

H2, 2 (p) = \[2{G^)f + (G 2 (2p)) 2 ], (51) 

H 2 M = Yl G{qf[2G{q) + G{q - 2p)\. (52) 

1 

The Fourier transformations are performed as FFT on one coordinate direction af- 
ter another and the whole 2 loop computation again has computational complexity 
DL A In L only, see [7] for more details. 
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With these expressions we can write (omitting the remainders. . . . + 0(</q)) 

9l Jisip's) ~ J2s{Ps) 



Y s = l + 



P's ~ PI 



(53) 



Am ( 2 = wig — m\ 



and 



9o T , 9o j tt ,9o 
— J u + jJuHu + - 



(54) 



Pt ~ P 2 t 



9s = 



9o -gl^H ls (p s ) +gl 



3H 2 ,l s {ps) + -H 2 ,2s{Ps) + 



( v , 1 J2s{p' s ) - -hsiPs 

2 h 2MPs) + 3- 



P's - PI 



(55) 



In order to obtain g s as a function of go and z t we have to combine now the last 
two equations to eliminate ml on the right hand sides. To the order considered 
and using 

—2 = - H lt{0), Jls , 2 = -2H 2 ,3s(p) 



dm\ dm\ 



we arrive at 



Am 2 = q x (t, z t , L)g + q 2 (t, z t , L)g 2 



with 



Qi(t,Zt,L) = ~-Ju(m 2 t ), 



q 2 (t,z t ,L) = - 



J2t(m 2 ,p t ) - ™£ + ^ t 2 (J 2t (m 2 : p' t ) - J 2t (m 2 ,p t )) 
Pt Pt 



and then at 



with 



9 = 9o+ Pi(s, t, z t , L)g 2 + p 2 (s, t, z t , L)gl 



(56) 
(57) 

(58) 
(59) 

(60) 



Pi(s,t,^,L) = --H ls (m 2 ,p s ), 

3 3 
p 2 (s,t,z t ,L) = 3H 2>ls (m 2 ,p s ) + -H 2>2s (m 2 ,p s ) + -H 2i3s (m 2 ,p s 



1 - 



J 



it 



J 



Is 



1 J 2s (m 2 ,p s >) - J 2s (m 2 ,p s ) 
3 



Pi ~ P 2 S 



(61) 
+ 

(62) 



In these formulae m t on the right hand sides is given by z t /L, of course, and m\ in 
the arguments of J ls , iJ ls , . . . refer to the mass value used here in the propagators 
that enter. 
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B.2 Couplings based on partition function ratios 

We here work out the coefficients appearing in (I3TI) . The leading order is trivially 
given by 

f/{z s ,L/a) = ~(£- E ) HP 2 + m 2 s ). (63) 

The first correction receives contributions from both the interaction term and from 
eliminating m 2 , for the renormalized mass m 2 s by (j54|) and we find after simple steps 

ff( Zs ,L/a) = ^-[G s (0)-G s (0)} 2 (64) 
with the mass m s here in both propagators. Numerical values are given in Table 



L 


Jo 


ft' A x 10 3 


f a,p 

Jo 


ft* x 10 3 


8 


0.26076127 


0.35227920 


0.30274830 


2.4978979 


10 


0.24731850 


0.33411586 


0.29837653 


2.4880152 


12 


0.24059900 


0.32490117 


0.29621943 


2.4832198 


16 


0.23423877 


0.31612031 


0.29419783 


2.4787956 


22 


0.23052145 


0.31096912 


0.29302625 


2.4762694 


32 


0.22836421 


0.30797464 


0.29234996 


2.4748253 


64 


0.22693150 


0.30598406 


0.29190233 


2.4738756 



Table 6: Perturbative coefficients fo' A , fi ,A for z a = 2. 



and asymptotic large L Symanzik expansions for these cases are 

/ a ' A = 0.226457 + 1.941L- 2 + 0(L~ 4 ), (65) 

ff A x 10 3 = 0.305324 + 2.698L" 2 + 0(L~ 4 ), (66) 

- / a,p = 0.291754 + 0.6049L- 2 + 0(L~ 4 ), (67) 

ft* x 10 3 = 2.47356 + 1.278L~ 2 + 0(L~ 4 ). (68) 
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